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’ DESIGN OF A MACH 8.0 AXISYMMETRIC NOZZLE FOR 

A HYPERSONIC TEST FACILITY 

SUMMARY 

An axisymmefric nozzle has been designed to produce test section 
flow at a Mach number of 8.0 for use in a hypersonic test facility 
at MSFC. Nominal design conditions used to calculate the viscous 
correction to the wall contour were selected from the results of a 
parametric boundary layer investigation so that the widest possible 
range of satisfactory operating conditions could be obtained. Coordi- 
nates for the nozzle are presented in tabular form suitable for design 
and manufacturing. The basic analysis techniques have been used to 
generate results to compare with experimental data from a facility at 
Langley Research Center. The agreement was reasonably good. 


INTRODUCTION 

An axisymmetric nozzle for a hypersonic test facility at MSFC has 
been designed to produce a test section Mach number of 8.0 when opera- 
ted at nominal design conditions. The nominal conditions were selected 
from the results of a parametric investigation of the effects of stag- 
nation pressure, stagnation temperature, and wall temperature upon the 
boundary layer correction to the inviscid nozzle coordinates. , 

These coordinates were obtained from a method-of-character is tics 
computer program devised at MSFC. This program, starting from an 
initial value line obtained* from a transonic solution and using a 
prescribed nozzle centerline velocity distribution as a boundary con- 
dition, computes the entire nozzle shape downstream of the throat. 

The transonic solution that was used is valid only for a circular arc 
approaching the throat from the subsonic side. Thus, the contour of 
a portion of the subsonic side is shaped by the transonic solution 
requirements. The remainder of the subsonic side was shaped to mini- 
mize the possibility of thermal stratification and to maintain com- 
patibility with the requirements of lower Mach number nozzles. 

A turbulent boundary layer program using an integral solution 
technique was used to obtain the viscous correction to the inviscid 
nozzle coordinates. The adequacy of the combination of these two 
programs has been demonstrated previously for lower Mach number nozzles. 
Included herein is a comparison of calculated boundary layer results with 
experimental data from the Langley Research Center Mach 6 hypersonic 
tunnel. The agreement between the two types of data is satisfactory. 
Since there should be virtually no Mach number effects on the accuracy 
of the inviscid solution, it is anticipated that the nozzle will pro- 
duce a satisfactory flow field for test purposes. 



ANALYSIS AND RESULTS 


The analysis that was applied to compute the supersonic side 
inviscid coordinates of the nozzle was the axisymmetr ic method of 
characteristics. In formulating a computer program to apply this 
method, which is a point-by-point solution of the entire supersonic 
portion of the flow field, initial conditions, centerline boundary 
conditions, and final Mach line boundary conditions must be specified. 
These conditions, as illustlated in sketch A, were specified in the 
following manner for the computer program used to generate the 
numerical results. 

The location of the initial value line and the numerical value 
of the flow variables along the line were obtained from the transonic 
solution in reference 1. Briefly, this solution is obtained by 
using a power series expansion in the equation of motion and solving 
for the unknown coefficients in terms of the radius of curvature of 
the nozzle wall, which is assumed to be a circular arc. This solution 
should be most accurate for nozzles which have a large radius of cur- 
vature compared with the throat radius. For this reason, a value of 

p„/r. = 10.0 was selected. 
a Z 

The centerline velocity distribution must be specified from 
point 1 (sketch A), whose location and velocity value are specified 
by the transonic solution, to point 2, which has a somewhat arbitrary 
location, but has the specified design velocity. This velocity dis- 
tribution may be specified in the program by a polynomial equation 
which has coefficients defined by input data or it may be defined 
point by point as input data. For this study, the latter procedure 
was selected since it offered two significant advantages in the 
hypersonic case: (1) The step size along the centerline can be 

variable in order to achieve maximum accuracy of the solution. 

(2) A velocity distribution can be obtained which restricts the 
maximum wall angle to any desired value. This technique can be used 
to obtain small wall angles to minimize start time. The major pro- 
blem with specifying the centerline velocity distribution in this 
manner is the additional time and effort it requires. 

The centerline velocity distribution used for this design study 
was obtained by computing the flow field in a nozzle which consisted 
of a circular arc downstream of the throat followed by a ten-degree 
conical nozzle tangent to the circular arc at the point of intersec- 
tion. This shape, shown in sketch B, was selected after a flow field 
analysis of several contours with different shapes leading up to the 
conical section. A value of 20 r t was selected for the circular arc 
radius of curvature; the flow field studies proved that this was 
sufficient to prevent the formation of a shock from the arc-cone 
juncture. 
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iONGI TWIN A L STATION - 





Calculations performed in the nozzle design deck require that the 
derivative of the centerline velocity be known. This was obtained by 
numerically differentiating the data obtained from the conical nozzle 
flow field solution. The final derivative function shown in figure 1 
was a composite of this curve and a linear, constant slope section 
tangent to the numerical function and extending to a value of zero for 
the derivative. Selection of the point of tangenc.y was governed by the 
requirement that the integral of the derivative function must yield the 
design velocity at the longitudinal station where the derivative goes 
to zero. The velocity distribution resulting from this composite 
derivative function, which was used in the design program, is'presented 
in figure 2. ’ 
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The boundary layer growth analysis is essentially the same as that 
derived by Elliot, Bartz and Silver (ref. 2). Their assumption of a 
l/7th power velocity and temperature profiles in the boundary layer has 
been replaced by ah arbitrary power determined from a correlation" of 
experimental data. 

■' ' f 

. Basically, the analysis is a simultaneous solution of the integral 
axisymmetr ic , compressible momentum and energy equations for an ideal 
gas.. Standard definitions of the boundary layer velocity, temperature, 
and displacement thicknesses, along with the assumption of an arbitrary 
velocity and total temperature power profile, are used to derive addi- 
tional equations to complete the boundary layer growth solutions:. >.ii. 
Constant pressure, is assumed across the boundary layer. v 

... • The boundary layer power profiles (velocity and total temperature) 
are obtained from a correlation of experimental data as a function of 
momentum Reynolds number (ref. 4). Friction coefficients are taken 
from a correlation of adiabatic, incompressible flat-plate experimental 
data and corrected for compressibility and heat transfer effects by use 
of Eckert's reference temperature (refs. 5 and 6). Stanton number is 
obtained from the Von Karman incompressible equation. 

As part of the design study, a boundary layer analysis was per- 
formed for the Langley Mach, 6 hypersonic tunnel configuration. Experi- 
mental data obtained from the calibration of this facility are available 
in reference 3. Comparisons of the analytical and experimental dis- 
placement ■; total boundary layer (velocity) arid 'momentum thicknesses 
are presented in figures 3, 4, and 5. Agreement between measured and 
analytical displacement and total boundary layer thicknesses (figures 
3 and 4) is seen to be good. However, agreement between the two sets 
of momentum thicknesses (figure 5) is only fair in the downstream 
portion of the test section and is relatively poor in the upstream 
section. 
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. Since /the, tunnel wall boundary layer description 1 at a given axial 
location defends upon upstream conditions and since experimental data 
were not obtained for the converging -diverging section of the Mach 6 
tunnel, the exact cause of the momentum thickness discrepancy could not 1 
be ascertained. However, there are two possible explanations. First, 
it was noted that the analytical and experimental velocity and tempera- 
ture profiles through the boundary layer' did not match in the upstream 
portion of the test section. This..* was .due, in part at le.ast., to the ' 
rionuniform temperature in the' stagnation chamber' reported' in reference 
2. Second, in calculating the analytical momentum thickness,, the 
integral momentum equation was ‘solved using a flatjplate friction co- 
efficient; whereas, the experimental data were obtained from the physi- 
cal, definition of momentum thickness using measured boundary layer' " 
properties.:' { 

■ Although- these deficiencies in the momentum thickness comparison 
were observed, the overall comparison of the analytical and experimental 
data is good. Thus, it appears that use of the analytical displacement 
thickness to correct the inviscid’ coordinates will provide a ’satisfac- 
tory tunnel wall contour. ‘ ? ' *•_ '■ 

Figures 6 and 7 present the Mach 8 axisymmetric nozzle boundary 
layer displacement thickness and total thickness, respectively. "These 
calculations are; based on a nominal set of stagnation chamber values; 
which are noted' on the' figured. The nominal design values 'were selected 
after a parametric boundary layer study over the entire range of stag- 
nation variables. Nominal values were selected that minimize the effect 
of deleterious boundary layer variations over the most useful operating 
range of the facility. 

The coordinates; for , the entire ; curvilinear portion of _ the nozzle 
wall are given in table 1. All of the coordinates are given in inches 
and the' throat has been used as the reference station/for the longitu- 
dinal coordinates. These coordinates, which are the end result of the 
entire. study, are. the values that will be used for nozzle fabrication. 
The subsonic shape from the valve to the first ordinate listed in 
table 1 consists of cylindrical and conical tubes and are explicitly 
defined by appropriate design drawings. 


CONCLUDING REMARKS 


: Analysis techniques incorporated in currently operational computer 
programs have been used to derive the nozzle contour for the design of 
a hypiersonic test facility at MSFC. These analysis techniques have 
been verified thrbugh comparison with experimental data and calibration 
results from other nozzles designed by the same procedures. Therefore, 
it is anticipated that the current design will produce a flow field 
of satisfactory quality. 
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